Spin relaxation and decoherence of two-level systems 
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We revisit the concepts of spin relaxation and spin decoherence of two level (spin-1/2) systems. 
From two toy-models, we clarify two issues related to the spin relaxation and decoherence: 1) For 
an ensemble of two-level particles each subjected to a different environmental field, there exists 
an ensemble relaxation time T* which is fundamentally different from Ti. When the off-diagonal 
coupling of each particle is in a single mode with the same frequency but a random coupling strength, 
we show that T* is finite while the spin relaxation time of a single spin Ti and the usual ensemble 
decoherence time T 2 * are infinite. 2) For a two-level particle under only a random diagonal coupling, 
its relaxation time Ti shall be infinite but its decoherence time T2 is finite. 

PACS numbers: 73.23.-b, 72.70.+m, 71.70.Ej 



Spin relaxation and decoherence is of central impor- 
tance in quantum computation and spintronicsQ, Q- 
Two phenomenological quantities, Ti and T2 known as 
the spin relaxation time and spin decoherence time for a 
two-level (spin-1/2) system, were introduced to describe 
the spin dynamics under the influence of an environment 
or external field in early theory As illustrated in 
Fig. HJa), a spin in its upper level, denoted by |+), may 
jump to its lower energy level, denoted by |— ), due to 
spin-environment coupling. T\ is defined as the average 
transition time. A quantum spin can also be in a super- 
position state, say (|+) + \ —))/V2 initially. Without the 
environmental interaction, the system undergoes a con- 
stant precession with a frequency u> = E/h, where E is 
the energy difference between |+) and |— ). Under the 
environmental influence, the phases of the expansion co- 
efficients in the superposition state will become random, 
and their relative phases, which describe the precession 
angle a in Fig. H^b), will be completely undetermined 
(decoherence). The time for the superposition state to 
lose its coherence is called T 2 . Ti = T 2 is believed to be 
true when the spin is isotropic while T 2 < 1T\ holds in 
general 00- 

In recent years, it has been realized that another quan- 
tity, , is important for a large number of experiments 
that measure the decoherence of an ensemble of spins 
rather than a single spin0,0|. Illustrated in Fig. Efc), m 
addition to the uncertainty in the precession angle of a 
single spin at a given time, there will be an uncertainty 
in the relative precession angles among different spins 
due to the spatial inhomogeneity. This additional deco- 
herence source leads normally to T 2 * < T 2 , within which 
T 2 * may deviate notably from T 2 . For example, T 2 * °f 
electron spins in quantum dots was shown to be several 



orders smaller than T 2 [ 
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FIG. 1: (Color) An illustration of Ti, T 2 and T 2 * . (a) Ti is the 
typical time for a spin jumps from the upper level (up spin) to 
the lower level (down spin). Inset: The polar angle j3 and the 
azimuthal angle a of a spin along unit direction n. (b) The 
precession of spin in xy-plane. At t — 0, the spin points to 
x-direction. Due to its irregular precession, the spin direction 
in a later time may be in a range of directions denoted by 
grey (blue) arrows. The spin loses its direction completely at 
T2, called the decoherence time, (c) For an ensemble of spins, 
all spins point to the x-direction at the beginning. As time 
goes on, different spins will precess to different directions due 
to the fluctuating field is different at different places. After 



Our current understanding 3] of T\ , T 2 and T 2 * s 
more accurate and much deeper than that of the early 
years^. Its applications to various systems, especially 
those of nanostructures are of current focus. However, 
there are still many vague conceptions or even miscon- 
ceptions about the description, origins and mechanisms 
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of these quantities. In this letter, we revisit the con- 
cepts of the spin relaxation and decoherence of two-level 
systems. We shall argue that, besides an ensemble deco- 
herence time T 2 , an ensemble relaxation time T* is also 
necessary to interpret experimental results. We show, 
through a toy-model calculation, that T* distinguishes 
itself from T\. In the case of our study, the total po- 
larization of the system decreases with time while each 
spin undergoes a well-defined periodic oscillation. Thus, 
there is a finite pseudorclaxation time T* for the whole 
system, but an infinite T\ . We investigate also the effects 
of a random diagonal coupling on spin relaxation and de- 
coherence. We show that both slow and fast fluctuation 
are harmful to T 2 while its effect on T\ must go through 
the off-diagonal coupling. 

Formulation of the dynamics of a two-level system. To 
mimic the motion of a two-level (spin- 1/2) system cou- 
pled to an environment, we assume that the system may 
jump from one state into the other through the coupling, 
or the system in a superposition state may obtain random 
phases to its expansion coefficients. The Hamiltonian of 
the system can be properly described by a 2 x 2 matrix, 



H=^a z + K(t) ■ a = 



± + K z K x -iK v 
K x + iK y -§-K, 



(1) 



er = cr x x + a y y + a z z is the Pauli operator. The 
environment-system interaction is approximated by a 
randomly fluctuating field K{t) — K x x + K y y + K z z, 
which could originate from the interactions of the spin 
with phonons, photons or with other surrounding spins. 
In general, the field is correlated both in time and in 
space. Its covariance measures the strength of the cou- 
pling which should also depend on the temperature be- 
sides of other factors. How to obtain the field from a 
microscopic Hamiltonian of both system and the envi- 
ronment is an interesting problem in both physics and 
chemistry [loj. For example, it has been shown [lH that 
electron-phonon, spin-orbit, and/or hyperfme interac- 
tions for electron spins in a quantum dot are equivalent 
to off-diagonal couplings. We will see that the diago- 
nal coupling K z (t) and the off-diagonal coupling K x ^(t) 
play different roles in determining T\ and Ti respectively. 
Without losing the generality, the energies of the upper 
and the lower levels (|±)) are assumed to be ±E/2, re- 
spectively, in the absence of the fluctuating field. The 
energy may be due to the Zeeman interaction for a spin 
in a constant magnetic field. 

Express the general wave function of our two-level sys- 
tem as ip(t) = a (t)\+) + b(t)\ — )i the equations for the ex- 
pansion coefficients a(t) and b(t), derived from the time- 
dependent Schrodinger equation, are (set h = 1) 

= [E/2 + K z ]a + [K x - iK y ]b, 
i^ = -[E/2 + K z }b+[K x + iK y ]a. (2) 



The issues are to find out the spin relaxation and deco- 
herence time for a given fluctuating field. Spin-relaxation 
time T\ can be obtained by examining how (|a(i)| 2 ) 
changes with time under an initial condition, where (. . .) 
denotes the average over the system state. T\ is the char- 
acteristic decay time extracted from (|a(i)| 2 ). 

To evaluate T 2 or T 2 1 if is useful to notice that a spin 
state along a particular direction n, specified by the polar 
and azimuthal angles (3 and a [inset of Fig. 1(a)], can also 
be written as ip(t) — exp(i(j))[cos((3/2) exp(ia/2)|+) + 
sin(/3/2) exp(— ia/2)\— )]. If one can find a{t) by solving 
Eq. (|5J) , then T 2 for a single spin or T 2 * f°r an ensemble of 
spins corresponds to the time at which the deviation of 
a equals it, i.e., (Aa 2 (T 2 )) = (a(T 2 ) 2 ) - (a(T 2 )) 2 = tt 2 . 

Ensemble spin relaxation time T* . In order to demon- 
strate the necessity of introducing an ensemble spin re- 
laxation time T* , we consider an ensemble of nonin- 
teracting spin-1/2 particles each of which is described 
by a Hamiltonian similar to Q with K z = 0, and 
K x + iK y = Ke lut .Then Eq.0 can be solved analyti- 
cally with solution 



a k {t) 
b k (t) 



e^ u M*[a k (0) cos(ffi) 
.A^a fe (0)-2^ fe 6 fe (0) . 
2TL sin(fii)], 

e i(w/2)t [6 fe (0) cos(ffi) 
.ALub k (0) + 2K k a k (0) . 



2Q 



sin(fit)], (3) 



whereAw = u — E and 4fi 2 = 4i^ 2 



Slu 2 



To make following analysis simple, we shall first as- 
sume uj = E for all spins, and K — K k > a random 
real constant for fcth spin. Thus, each spin is nothing 
but the famous Rabi problem at resonance. For the kth 



spin, Eq. © gives a k {t) = e 
ib k (0)sm(K k t)} and b k (t) 



= p -*(E/2)t 
J{E/2)t 



a k (0)cos(K k t) - 
[b k (0)cos(K k t) - 



JOfc(O) sin(K k t)]. If one measures the expectation value 



of S z of the fcth spin defined as (S k 



\b k (t)\ 2 \ (His set to I), one finds (S k<z ) = [\a k (0)\ 2 - 
\b k (0) | 2 ]/2 cos(2K k t) + %[a k (0)b* k (0)] sm(2K k t) , which os- 
cillates periodically in time instead of decay. One then 
should conclude T\ = 00 as it should be. However, if 
we measure S z of the whole system, wc need to include 
contributions from all spins. Since different spins have 
different K values, they will oscillate with different peri- 
ods, and the spins initially in phase will be out of phase 
during the evolution, as is illustrated in Fig. 2. Thus, 
even all spins are fully polarized along z direction ini- 
tially, S z will decrease with time. If K has a distribution 
function of P{K) — (2K/a)e~ K where a measures 
the width of the distribution, by using Eq. 10) , we have 

[(Sk)] = (l/N)Ek(Sk,z) = 1~ V^te~° t2 erfi(^dt), 
where erfi(x) is the imaginary error function, and [• • • ] 
is an average over all spins. This time dependence is 
plotted in the inset of Fig. 2. Thus, one finds the polar- 
ization of the system relaxes after a characteristic time 
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Ty = 1/ y/a, although there is no relaxation for individ- 
ual spin. It is also interesting to evaluate T 2 and T 2 of 
this toy model. Set afe(O) = &fe(0) = 1/V2 for all k and 
compare cik{t) and &&(£) with their expressions in terms 
of the polar and azimuthal angles mentioned in the early 
paragraphs, we find afc(i) = Et. All spins precess at the 
same rate and in phase. Thus, both T 2 and T 2 * are infi- 
nite. This is a peculiar case in which each spin follows 
a deterministic motion, and there is no relaxation and 
no decoherence. But if one measures the polarization of 
the whole system, the system appear to have a relaxation 
time T* that is fundamentally different from T\. 

In fact, one may also assume that fcth spin has lu^ 
and different level spacing Ek- Since E^, Uk, and Kk 
are constants for a given spin, a&(t) and bk{t) are well- 
defined periodic functions according to Eq. @ . Thus its 
precession angle and expectation value of Sk, z , which can 
be obtained after some algebras 



MQ)| 2 -k(Q)|- 

2ft 2 



M/4 + #2cos(2fit)], 



which oscillates periodically in time. Therefore, both 
T\ and T 2 are infinite. However, if one measures S z or 
S x — iS y (which is a measure of precession) of the whole 
system, we need to sum over all spins. This corresponds 
to a sum of many periodic functions with different peri- 
ods. Similar to the case of random Kk discussed early, 
it leads to a decay function, and T* and T 2 *are finite be- 
cause different spins precess with different speeds in this 
case. 
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FIG. 2: The time dependence of occupation probability of 
different spin in its upper level. All spins are initially in their 
upper levels. Due to the random coupling strength, their Rabi 
frequencies are also random. Inset: The time dependence of 
[{Sk)]. 

Away from the resonance spin initially at the 

upper level will not evacuate completely from the level. 
The occupation probability for a spin in the upper level 
is given by the well-known Rabi solution 



where 4ft 2 = AK 2 + (w - E) 2 , and A(u) = AK 2 /[AK 2 + 
(u> — E) 2 ]. Amplitude A at the resonance takes the max- 
imal value 1. A decreases fast when ui deviates from E. 
Thus, the mode of the fluctuating field near the resonance 
is expect to relax the spin more effectively. We shall see 
that the diagonal coupling behaves differently. 

Effects of diagonal coupling on T± and T%. We would 
like to discuss how a diagonal coupling affects the spin 
decoherence and the spin relaxation. Consider a case in 
which K z in Eq. Q is a randomly fluctuating field and 
K x = K y = 0, this is an exactly solvable problem because 
there is no coupling between different levels, and K z 
causes only the energy levels to move up and down ran- 
domly. The consequence is the accumulation of a dynam- 
ical random phase for each level. The solutions for the 
expansion coefficients are a(t) = a(0)exp[— i J Q (E/2 + 

K z {t'))dt'\ and b(t) = 6(0) exp[i J*(E/2 + K z (t'))dt'], 
where a and b are separated because of the absence of 
off-diagonal matrix elements. Thus, a spin in either up- 
per or lower level will not relax, and T\ is infinite. How- 
ever, for a spin in superposition state, each of the two 
expansion coefficients will obtain a random phase due to 
random field K z , leading to a random precession angle 
a(t) = Et + 2 J* K z (t')dt'. Thus, the expectation of the 
deviation of the precession angle is 



(Aa 2 (t))=4 (ftidt 2 (X*(ii)X*(i a )). (5) 
Jo 

For a stochastic field with autocorrelation, 

(K z (tx)K z (t2)) = ^exp(- 1 -^^), (6a) 



\a(t)\ 2 = 1 - A{uj) sin 2 (fit), 



(4) 



where r c and A are the correlation time and strength, we 
have 

(Aa 2 (t)) — 4A[f + r c (e~^ — I)]. (6b) 

Thus, the spin-decoherence time T 2 is given by 

T 2 + r c [exp(-T 2 /r c ) - 1] = tt 2 /(4A). (7) 

In the limit of r c — > 0, the fluctuating field is a white noise 
{K z (t\)K z (t2)) — A6(ti — ^2)7 and the decoherence time 
is simply T 2 = 7r 2 /(4A), inversely proportional to the 
correlation strengthen of fluctuating field and indepen- 
dent from the level spacing. Therefore, if the amplitude 
of the fluctuating field increases n times, T 2 decreases by 
n 2 times. 

Solution exhibits three regimes. In the white noise 
regime where r c A <C I, the fluctuating field behaves 
like an independent force on the spin precession. The 

spin precession angle undergoes a random walk around 

2 

its mean rate da/dt = E, and T2 fts j^, insensitive 
to t c . In the long correlation regime where the spin 
experience a random persistent force, the deviation of 
the precession angle increases quadratically with time, 
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and the motion of the spin behaves like a biased ran- 
dom walk. Correspondingly, decoherence time T 2 is 
T 2 « 7Ta/ r c /(2A) <C r c . In the intermediate regime of 
r c A ~ 1, the T 2 , r c and A -1 are in the same order. Un- 
like the off-diagonal coupling that affects the spin relax- 
ation most when it is at the resonance, all fluctuations 
of the diagonal coupling contribute to the deviation of 
the precession angle as it is shown in Eq. Thus, 
no matter whether there are fast (white noise) or slow 
(long correlation time) fluctuations, K z will lead to an 
effective decoherence. Of course, it does not mean that 
all fluctuations will play an equal role. 

Our analysis so far does not rely on the sources of the 
fluctuating field. Thus, we expect our results to be ap- 
plicable to all two level systems as long as they can be 
described by a Hamiltonian similar to Q). Although T* 
in our toy-model is more a definition in a thought exper- 
iment than a measure of a true spin relaxation, it makes 
the point that the relaxation time for a single spin and 
for an ensemble of spins could be totally different. We 
present an extreme case where Tj" is finite while T%, T 2 
and T 2 are infinite. In many realistic two-level systems, 
one may find another extreme where the distinction be- 
tween Ti and T* is not necessary but the existence of 
realistic systems between the two extremes so that one 
has to use two relaxation time to separately describe the 
relaxation for a single spin and for an ensemble of spins, 
just as To and T 2 * for spin decoherence, cannot be ruled 

outHlli. 

It may be proper to explain the processes that con- 
tribute to Ti (I?) and T 2 (T|). T x and T* are caused by 
the transitions of a system from one quasi energy level to 
another. Such a transition is usually irreversible, accom- 
panying a complete change of the wave function. This 
microscopic view of Ti, T* is not in contradiction to 
the thermodynamic view that Tiis the time of thermal 
equilibration for spin population 3} because the time for 
a system to reach its thermal equilibrium normally re- 
quires many jumps between the upper and lower levels. 
Of course, the numbers in the two definitions may be 
different, but they should be in the same order. This is 
similar to the definition of the electron relaxation time 
scattered by impurities, which is the characteristic time 
for the electron distribution to return to its equilibrium 
one; on the other hand, it can be calculated from the 
Fermi golden rule at the quantum mechanical level. T 2 
and T 2 come from the loss of coherence of a state which 
could be due to a transition for Ti (T*), but it could 
also due to the accumulation of random phases to the 
components of a superposition state without a transi- 
tion between two levels, as we showed in our toy-models. 
Thus, in general, T 2 and T 2 * must be smaller than Ti 
and T*. In our extreme case, the T£ is not caused by 
true irreversible transitions, but by coherent transitions 
(Rabi oscillations). In this sense, T* in our example is a 
pseudorelaxation time, and it should not be surprising to 



see our T* deviating from T 2 < T*. There are miscon- 
ceptions in the literatures about the difference of T 2 and 
T 2 * . It was claimed that T 2 is referred to as the reversible 
dephasing 3], while irreversible dephasing was attributed 
to T 2 . This claim is questionable because the nature of 
randomness is irreversible. Thus both T 2 and T 2 are due 
to the irreversible loss of coherence. If the irreversibility 
in the claim is referred to the transition from one energy 
level to another, it is again not true because both T 2 and 
T 2 * can be caused by the accumulation of random phases 
on the expansion coefficients, as revealed in our second 
toy-model. 

Due to the different natures of Ti and T 2 , the diago- 
nal and off-diagonal coupling play quite different roles on 
the spin relaxation and the spin decoherence. First, with 
only diagonal coupling, there is no spin relaxation. But 
there is a spin decoherence. However, a diagonal coupling 
can also influence the spin relaxation if an off-diagonal 
coupling is present. For an off-diagonal coupling, those 
fluctuation modes near the Rabi resonance affects spin 
relaxation more. Diagonal coupling could affect the re- 
laxation either through the shifting of resonance point or 
through the modification of quantum interference by the 
phase change of wave function. 

It may also be useful to comment on the relationships 
between T 2 and T 2 . For a non-interacting spin system 
subject to the same fluctuating field, we should have 
T 2 < T 2 because, in addition to the causes of spin deco- 
herence on individual spin, the spatial inhomogeneity can 
contribute to decoherence of an ensemble of spins. The 
importance of the inhomogeneity on the decoherence was 
known in many fields and also in early work in semicon- 
ductor community, [jl 0] in which the interference effect 
due to the inhomogeneity of k space or of real space can 
cause decay. For an interacting spin system, it is more 
complicated. An interaction can be an extra source of a 
fluctuating field. On the other hand, the interaction may 
glue different spins together. The bondings among spins 
could be so strong that all spins have to move coher- 
ently. In this case, one can expect that the dynamics of 
the system slows down. Also, a system with a single spin 
and a system with many spins are totally different, and 
there is no comparison. For example, the spin relaxation 
and spin decoherence of a Co atom is completely differ- 
ent from that of a bulk Co magnet where all spins are 
aligned in the same direction. And it is not surprising to 
have spin decoherence T 2 * of the Co-magnet to be much 
longer than T 2 of one Co-atom. Another example may 
be superconductor caused by the effective e-e attraction. 
It is known that coherence time of electrons in a super- 
conducting state could be much longer than that of an 
electron in a metal. 

In conclusion, we demonstrate the necessity of intro- 
ducing an ensemble spin relaxation time T* that could 
be orders of magnitude different from Ti . We have also 
shown that all fluctuations should be relevant to the spin 
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decoherence. There is no too-slow or too-fast fluctuation 
for T 2 and T 2 *. 

This work is supported by UGC, Hong Kong, through 
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